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Abstract 

Studies of leukocyte telomere length (LTL) and adiposity have produced conflicting results, and the relationship between 
body mass index (BMI) and telomere length throughout life remains unclear. We therefore tested association of adult LTL 
measured in 5,598 participants with: i) childhood growth measures (BMI and age at adiposity rebound (AR)); ii) change in 
BMI from childhood to adulthood and iii) adult BMI, waist-to-hip ratio (WHR), body adiposity index (BAI). Childhood BMI at 
AR was positively associated with LTL at 31 years in women {P = 0.041). Adult BMI and WHR in both men {P = 0.025 and 
P = 0.049, respectively) and women (P = 0.029 and P = 0.008, respectively), and BAI in women (P = 0.021) were inversely 
associated with LTL at 31 years. An increase in standardised BMI between early childhood and adulthood was associated 
with shorter adult LTL in women (P = 0.008). We show that LTL is inversely associated with multiple measures of adiposity in 
both men and women. Additionally, BMI increase in women from childhood to adulthood is associated with shorter 
telomeres at age 31, potentially indicating accelerated biological ageing. 
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Introduction 

Telomeres are protective DNA-protein structures that cap the 
ends of Knear chromosomes, extending to 10-15 kb in humans. In 
proliferative tissues such as leukocytes, telomeres shorten with each 
cell division - a process believed to be accelerated by oxidative 
stress and inflammation [1—3]. Mean leukocyte telomere length 
(LTL) in adulthood is influenced by genetic factors [4], and 
multiple environmental exposures including chronic psychological 
stress [5], childhood adversities [6,7] and lifestyle factors such as 
smoking [8] . LTL is positively correlated with human lifespan [9] , 
while shorter LTL is associated with several age-related conditions 



including type 2 diabetes (T2D), heart disease and some cancers 
[10-13]. LTL has thus been proposed as a biomarker of biological 
ageing, though further longitudinal studies are required to confirm 
this hypothesis [14,15]. 

Inverse associations between telomere length and body mass 
index (BMI), waist-to-hip ratio (WHR), waist circumference and 
visceral fat in adulthood have been reported [16-20]. However, 
other studies have indicated a lack of association with these 
adiposity-related phenotypes (reviewed in [21]). Furthermore, the 
relationship between telomere length and body adiposity index 
(BAI, calculated as (hip circumference (cm))/((height(m))' "')— 18) 
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has not been investigated, despite assertions that this anthropo- 
metric measure may be a better index of body fat than BMI [22] . 
One explanation for the conflicting results for associations between 

LTL and measures of adiposity could be sample diversity both 
within and between studies, in particular the inclusion of wide age- 
ranges and uneven sex distributions. Indeed, one study reports a 
stronger inverse association between BMI and LTL in young (< 

30 years) compared to older (>60 years) women [16], suggesting 
that the relationship between obesity and telomere length may 
change throughout life. In recognition of these inconsistencies in 
the current literature, a recent review and meta-analysis of 
previous studies of telomere length and BMI called for further 
large-scale epidemiological studies in different age groups and 
sexes, ideally including additional measures of adiposity [23]. Our 
objective was to address this gap in research. 

With regard to childhood obesity, we previously carried out a 
case-control study in a cohort aged 2-17 years. This analysis 
revealed that severely obese children have shorter white blood cell 
telomeres in childhood than their normal weight peers [24]. 
However, to our knowledge, no workers have investigated whether 
childhood BMI, or the increase in BMI from childhood to 
adulthood, is associated with telomere length in adulthood. Given 
the rising prevalence of both child and adult obesity worldwide, 
and the potential increased risk of age-related conditions such as 
T2D and cardiovascular disease that is associated with excess or 
abnormal fat accumulation, these are key issues to address [25,26]. 
The biological mechanisms that may underlie associations 
between shorter LTL and increased BMI are currently unknown, 
although obesity is recognised as a state of both increased oxidative 
stress and low-grade inflammation: both processes that are 
believed to accelerate leukocyte telomere attrition [1-3]. 

In the present study, we aimed to investigate the relationship 
between telomere length and adiposity using a longitudinal design, 
making use of prospectively gathered data from the Northern 
Finland Birth Cohort 1966 (NFBC 1966). W(; us(xi childhood 
growth data to model age and BMI at adiposity rebound (AR), i.e. 
the point at which BMI reaches its nadir in early childhood 
(~5 years). Early AR has been shown to predict increased 
adiposity in adolescence [27] and may also be related to poorer 
metabolic health in adulthood [28]. Using anthropometric and 
telomere measurement data obtained for 5,598 participants aged 

31 (48% male), our aims were threefold: 

i) To investigate the association between age/BMI at AR in 
childhood and LTL at age 31; 

ii) To investigate the association between longitudinal change in 
BMI from childhood (~5 years) to age 31 and LTL; 

iii) To test the association between WHR, BMI, BAI and LTL at 
age 31. 

Materials and Methods 

Study Samples 

Pregnant women hving in the provinces of Oulu and Lapland 
with expected delivery dates in 1966 were invited to participate in 
the NFBC 1966, a prospective follow-up study from birth into 
adulthood [29] (Figure 1). A total of 12,058 live-born children 
(6,169 boys and 5,889 girls) were enrolled in the study; 96.3% of 
all births in the region. Detailed information on participants at 
birth, including mother's parity and parents' socio-economic status 
(SES), was collected using questionnaires at hospitals or by trained 
nurses. Childhood growth measurements were obtained by trained 
nurses from birth to 16 years of age. At age 31 years, all those alive 



and with known address were traced (N = 1 1,541) and sent postal 
questionnaires, from which information including participant's 
children, SES, smoking status and age at menarche was obtained. 
For this study, smoking status at 31 years was defined as either 
non/light (<10 cigarettes per day) or heavy (>10 cigarettes per 
day), since there was no significant difference between the mean 
telomere length of non-smokers and light smokers (data not 
shown). A subset of 8,463 living in the original target area or in the 
Helsinki area were also invited to a clinical examination, in which 
6,007 consenting participants took part [30]. Blood samples were 
taken at that time, from which 5,753 DNA samples were 
successfully extracted. LTL measurements were obtained for 
5,620 of these individuals, of which 22 extreme values were 
excluded (see below), leaving 5,598 for use in analyses. 

Ethics Statement 

Informed written consent for the use of the data and DNA was 
obtained from all subjects, and approval granted by the Ethics 
Committee of the Northern Ostrobothnia Hospital District in 
Oulu (Finland). 

Childhood Anthropometric Measurements 

For each individual, predicted age and BMI at AR were 
calculated using sex-specific polynomial mixed effects models, 
fitted on BMI measurements collected longitudinally from 
18 months to 13 years of age, with the mid-point at 7.25 years. 
A natural logarithmic transformation of BMI was used to reduce 
skewness. A detailed description of model selection and fit is given 
elsewhere [31]. The main model used was as follows: In (BMI (kg/ 
m^)) = P„+ Pi Age + P2 Ageh P3 Age'^+ P4 Sex + P5 Age * Sex + 
Pe Age^ * Sex + U()+ ui (Age) + s. 

The Po, Pi, P2, Ps, P4, Ps and Pg are the fixed effects, uq and ui 
are the individual level random effects and E is the residual error. 
Age at AR was defined as the age when BMI measurement 
reached its minimum, between 2.5 and 8.5 years, with cut-off 
points chosen from descriptive analysis of growth curves. Each 
individual's data were weighted by the number of measurements 
within the age window to account for uncertainty in the derived 
parameters. Individuals with fewer than three measurements were 
excluded from the analysis, as were multiple births. 

Adulthood Anthropometric Measurements 

BMI at 31 years was calculated from weight and height 
measurements using the formula weight/height^ (kg/m^). Individ- 
uals with a BMK18.5 were defined as underweight, those with a 
BMI>18.5<25.0 were defined as normal weight, those with a 
BMI>25.0 were defined as overweight and those with a BMI> 
30.0 were defined as obese [32]. 

WHR was calculated from waist and hip measurements taken at 
age 31 years. BAI at 31 years was calculated as described by 
Bergman et al [25], as (hip circumference (cm))/((height(m))''^) — 
18) (Table 1). 

Longitudinal Change in BMI 

Firsdy, values for zscore BMI at AR in childhood and zscore 
BMI at 3 1 years were calculated separately in males and females, 
standardised to zero mean and unit variance. The change in zscore 
BMI was then calculated as the difference between these two 
measurements: ((BMI zscore at 31) — (BMI zscore at AR)). The 
resulting variable, BMI zscore change, provides a sex-specific 
measure of the standardised increase in BMI from childhood to 
adulthood for each participant. 
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Total live born infants in NFBC 1 966 
N=12,058 



Target population for clinical examination were 
those living in the original area or in Helsinki 
N=8,463 



Deceased or address 
not known 

N=517 



Did not participate 
in clinical 
examination 

N=2.430 



Participated in clinical examination and 
consent for data use 
N=6,033 (71.3%) 






\ 






No consent to use 
data 
N=26 





Data from clinicnl examination. N=6.007 



Individuals with 
missing (N=133) 

or excluded 
(N=22) telomere 
measurements 



DNA extracted 

N=5.753 



Smgletons with growth data 
and consent for data use 
N= 4.121 



Indi\iduals with 

telomere 
measurements 
N= 5.598 



Childhood BMI at Adiposit%- 
Rebound (AR) from growth 

modelling 

N= 3.921 



Individuals with 

less than 3 
measurements 
for growth 
mode nine 
N=200' 



luduiduals with 
missing adult BMI 
(N=92). or 
excluded due to 
extreme adult BMI 
or BMI z-score 
change values 
(N=102) 



Three datasets: 

1) Individuals with telomere measurements and 
age BMI at AR in childhood. N=3.903 

2) Indi\iduals with telomere measurements and adult 
BMI (at 31 years). N=5.404 

3) Individuals with telomere measurements and BMI 
z-score change (child-adult). N=3.683 



Individuals with 
excluded telomere 
measurements 
N=18 



Figure 1. Flow chart for inclusion of study participants. 

doi:10.1371/journal.pone.0099133.g001 
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Leukocyte Telomere Length (LTL) Measurements 

Mean relative LTL was measured in genomic DNA samples 
prepared from peripheral blood leukocytes taken at age 3 1 , using a 
multiplex quantitative real-time PGR method [33], with minor 
modifications as described previously [24]. Briefly, the multiplex 
qPCR method for measuring mean relative telomere length 
provides a "T/S ratio" for each DNA sample. This is a relative 
measure of the amplification of the telomeric DNA sequence (T) 
compared to that of a single copy gene (S) in each test sample, 
normalised using a common reference DNA sample. The mean 
values for the calibration curves, based on values obtained for 
serial dilutions of the reference sample spanning 5-50 ng and run 
in triplicate, were 0.94 (SD = 0.02) and 0.96 (SD = 0.01) for the T 
and S amplicons respectively. 

The overall mean coefficient of variation (C V) for T/S values of 
duplicate test samples on the same plate was 5%, and the mean 
inter-run CV for selected samples was 6.2%. 

Exclusions and Sample Sizes 

Individuals with extreme BMI values (BMI>50 kg/m^ or a 
BMI 7.-score change >3 standard deviations from the sex-stratified 
population means, N= 102) were excluded to reduce risk of 
misclassification due to measurement or data entry errors. 
Additionally, samples with extreme LTL values were excluded 
from the analysis (>3 standard deviations from the overall 
population mean, N = 22). This relatively conservative cut-off 
was chosen to minimise potential uncertainty about the validity of 
extreme telomere length values obtained using qPCR [34]. The 
final number of participants investigated in each part of the study, 
after exclusions and taking account of all available data, is shown 
in Figure 1. Li summary, of the N = 5,753 participants who 
participated in the clinical examination at age 31 and had DNA 
samples available, three datasets were selected for analysis. 1) 
individuals with telomere measurements at 3 1 years and childhood 
BMI measurements at AR, N = 3,903; 2) individuals with telomere 
measurements at 31 years and adult BMI measurements (at 31 
years), N = 5,404; and 3) individuals with telomere measurements 
at 31 years and BMI i-score change (age at AR to 31 years), 
N= 3,683. The samples used for analysis were representative of 
the NFBC1966 subjects for these parameters and there were no 
difierences observed between the complete dataset and incomplete 
datasets (P>0.05). 

Statistical Analyses 

Descriptive characteristics of the study samples were determined 
overall and for men and women separately (Table 1). Means with 
standard deviation (SD) were calculated for continuous variables, 
and percentages were calculated for categorical variables, as well 
as for individuals classed as underweight, normal weight, 
overweight and obese. Student's t-test for unpaired data was used 
to test for difierences between males and females, and a Chi- 
square test was used to test for a difference in percentage 
distribution in categorical data, with a Pjifl- <0.05 considered 
significant. Kolmogorov-Smirnov tests were performed to deter- 
mine whether distributions for continuous variables followed 
normal distributions, and if not then variables were transformed 
accordingly. 

Relative LTL measurements (T/S ratios) were log transformed 
to achieve normality for use in sex-stratified analyses. LTL was 
considered as the outcome variable, while age and BMI at AR; 
WHR, BMI, BAI in adulthood and the BMI zscore change were 
considered as predictor variables. Two models were tested: 1) an 
unadjusted model, testing for association between log-transformed 
LTL and adiposity-related measurements from childhood and 
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adulthood; 2) an adjusted model, testing for association between 
log-transformed LTL and adiposity-related measurements from 
childhood and adulthood, controlling for potential confounders 

known to be associated with either the outcome or predictor 
variables. Those selected were maternal parity and SES at birth 
(or paternal SES at birth if maternal SES was missing), and age at 
menarche (females only), SES at 31 years, adulthood smoking and 
number of children at 31 years. To control for possible batch 
effects in the telomere measurements, all analyses were also 
adjusted for qPCR plate. Estimated effects and 95% confidence 
intervals are reported as % change in LTL per unit change in 
predictor variable. The Benjamini and Hochberg method was 
used to address the issue of multiple testing [35], with an adjusted 
P value (Pcmciai) of <0.05 considered statistically significant. All 
statistical analyses were performed using SAS, version 9 (SAS 
Institute Inc., SAS Campus Drive, Gary, North Carolina 27513, 
USA). 

Results 

Cohort Characteristics 

Table 1 shows the descriptive statistics ioi measures of 
adiposity, LTL and potential confounding varialjlcs, (or the study 
population overall and stratified by sex. Consistent with other 
human population studies [36], we observed a positive association 

between telomere length and female gender, with men having a 
mean T/S ratio of 1.12 and women a value of 1.17 (P<0.001). 

Association between Childhood Growth and Adiposity 
Measures and LTL at Age 31 

BMI at AR was positively associated with LTL at 31 years in 
females in unadjusted analyses (P= 0.007), and after adjusting for 
potential confounding factors and correcting for multiple testing 
(■fcorraferf^ 0.041). There was no association between BMI at AR 
and LTL at 31 years in males (Table 2). There was no evidence 
for an association between age at AR and LTL at 31 years in 
either males or females, in unadjusted analyses (Model 1) and after 
adjusting for potential confounding factors (maternal parity, SES 
at birth, smoking, SES and having children at 31 years, and 
additionally, in females only, age at menarche) (Model 2). 

Association between Adulthood Adiposity-related 
Phenotypes and LTL at Age 31 

In the unadjusted Model 1, WHR was inversely associated with 
LTL at 31 years in both males {P= 0.014) and females (P<0.001). 
After adjustment for potential confounders (maternal parity and 
SES at birth, age at menarche (female participants), SES, smoking 
and number of children at 31 years) and correction for multiple 
testing (Model 2), the association remained significant in females 
[Pamnckd = 0.008) and males (P,„r„,cii!d - 0.049). An inverse 
association between BMI and LTL at 31 years was also observed 
in both males {P— 0.002) and females {P- 0.005), which remained 
significant in males (/"^o^^^j,^ = 0.025) and females (/"con-ecferf^ 0.029) 
after adjustment for potential confounders and correction for 
multiple testing (Model 2). Body fat percentage, estimated by BAl 
at age 31 was also inversely associated with LTL in females 
{P- 0.004) and remained significant after adjustments 
(^correcferf— 0.021). However, BAI was not associated with LTL in 
males (Table 2). 



Association between BMI z-score Change from 
Childhood-adulthood and LTL at 31 

BMI zscore change from AR in childhood (at around 5 years) to 
3 1 years was inversely associated with LTL at 3 1 years in females 
(P<0.001) but not males (P= 0.238) (Model 1). A one unit increase 
in BMI zscore change was associated with a 3.08% decrease in 
LTL in females (Model 1), which reduced slightly to 2.91% after 
adjustment for potential confounding factors {Pcomcted— 0.008). 

Discussion 

We have shown in a large population-based cohort that a 
greater increase in standardised BMI from early childhood to 
adulthood is associated with shorter mean relative LTL in women 
aged 3 1 . We have also shown that BMI and WHR at age 3 1 are 
inversely associated with adult LTL in both men and women in 
this study, and that higher BAI is inversely associated with LTL in 
women only. AH these associations remain significant after 
adjusting for multiple potential confounders (maternal parity and 
SES at birth, and age at menarche (females only), SES, smoking 
and number of children at 3 1 years) and correcting for multiple 
testing. Our results thus confirm and extend those reported in a 
recent meta-analysis of BMI and LTL [23]. 

Our finding that WHR and BMI are both associated with LTL 
in adult women is consistent with previous reports of an inverse 
relationship between measures of adiposity and concurrent 
telomere length [8,16,19,20,37]. In contrast to Nordljall et al 
[37], we found that BMI is associated with shorter telomeres in 
adult men as well as in women. In addition, we identified an 
inverse association between WHR and LTL in adult men, which is 
consistent with a previous report that WHR is a significant 
predictor of faster telomere shortening rate in a predominantiy 
male cohort (>80'M)) [17]. To our knowledge, our study is the first 
to show that LTL is also inversely associated with BAI in women. 
We did not see any evidence for an association between BAI and 
LTL in men, perhaps because BAI may be more strongly 
correlated with percentage of body fat in women than in men [38]. 

The gender dilference for our main finding - that a greater 
increase in standardised BMI from childhood to 31 years is 
associated with shorter adult LTL - may reflect the fact that a 
greater proportion of change in BMI in women compared to men 
over this period is characterised by increased body fat [39] . This 
association remains significant even after adjusting for number of 
children born to participants by age 3 1 years, suggesting that this 
finding is unrelated to pregnancy. 

We did not find any evidence for an association between age at 
adiposity rebound (AR) and LTL at age 31 in either sex, despite 
previous studies in this cohort and others showing that earlier age 
at AR predicts higher BMI and poorer metabolic health in 
adulthood [28,40]. Intriguingly, we identified 2l positive association 
between BMI at AR and adult telomere length in women, which 
remained significant after correcting for maternal parity and SES 
at birth, and age at menarche, SES, smoking status and number of 
children at 31 years. Furthermore, there is no evidence of 
association between BMI at AR and cardiometabolic profiles at 
age 31 in this cohort [40]. We speculate that these findings may 
reflect a beneficial efiect of having a childhood BMI within the 
normal range (as did the majority of participants in this population 
cohort) on adult telomere length and metabolic health. Further 
longitudinal studies of growth throughout childhood, weight gain 
in adulthood and telomere measurements taken at multiple time 
points are required to elucidate the mechanisms underlying these 
observations. 
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The major strengths of our study include the uniform age and 
ethnicity of the participants, the even sex ratio and the availabiUty 
of LTL data for a large number of samples. Although concerns 
have previously been expressed over the reliability of the qPCR 
method [34], the telomere measurements used in this study show 
the expected sex difference, and have previously been used in the 
identification of robust genetic associations with LTL [4]. 
Additionally, the use of growth modelling phenotypes allowed 
the investigation of the relationship between important measure- 
ments of childhood adiposity and LTL in adulthood. 

Due to the availability of DNA samples for only one time-point, 
it was not possible to study the association between longitudinal 
telomere length change and BMI trajectories over the same 
period. Thus, we were unable to draw any conclusions regarding 
causality for the relationship between BMI increase and telomere 
length. We were also unable to study the potential effects of 
interim changes in BMI between childhood and adulthood. 
Finally, it would have been interesting to further investigate the 
effects of underweight and obesity in childhood and adulthood on 
adult LTL, but we did not have sufficient numbers of such 
individuals available in this population cohort. 

Despite these limitations, our study adds to the growing 
evidence supporting the hypothesis that factors associated with 
obesity, such as increased systemic levels of oxidative stress and/or 
inflammation, may have similar effects on lifespan and cellular 
function to the ageing process [41]. Our findings warrant further 
investigation of the potential for using telomere length as a 
biomarker to assess the effectiveness of interventions to combat the 
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negative impact of obesity on healthy ageing. These include 
increased moderate physical activity, which has been associated 
with longer telomeres (reviewed in [42]), as have dietary and 
holistic weight loss interventions [43,44]. 

Conclusion 

We show that LTL is inversely associated with multiple 
measures of adiposity in both men and women, and that BMI 
increase in women from childhood to adulthood is associated with 
shorter telomeres at age .31. Further longitudinal studies are now 
required to investigate the long-term effects of weight gain after 
childhood on LTL, healthy ageing and disease risk, and to 
investigate the gender differences observed. 
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